Complete Parameterizations of Diffusive Quantum Monitorings 



O 

(N 

(N 



Ph. 



> 

o 

m 

O 



X 



A. Chia and H. M. Wiseman 
Centre for Quantum Computation and Communication Technology (Australian Research Council); 
Centre for Quantum Dynamics, Griffith University, Brisbane, Queensland ^i-Zi, Australia 

(Dated: January 12, 2013) 

The master equation for the state of an open quantum system can be unraveUed into stochastic 
trajectories described by a stochastic master equation. Such stochastic differential equations can 
be interpreted as an update formula for the system state conditioned on results obtained from 
monitoring the bath. So far only one parameterization (mathematical representation) for arbitrary 
diffusive unravellings (quantum trajectories arising from monitorings with Gaussian white noise) of 
a system described by a master equation with L Lindblad terms has been found [H. M. Wiseman and 
A. C. Doherty, Phys. Rev. Lett. 94,070405 (2005)]. This parameterization, which we caU the U-rep, 
parameterizes diffusive unravellings by + 2L real numbers, arranged in a matrix U subject to three 
constraints. In this paper we investigate alternative parameterizations of diffusive measurements. 
We find, rather surprisingly, the description of diffusive unravellings can be unified by a single 
equation for a non-square complex matrix if one is willing to allow for some redundancy by lifting 
the number of real parameters necessary from + 2L to SL^ + L. We call this parameterization 
the M-rep. Both the M-rep and U-rep lack a physical picture of what the measurement should look 
like. We thus propose another parameterization, the B-rep, that details how the measurement is 
implemented in terms of beam-splitters, phase shifters, and homodyne detectors. Relations between 
the different representations are derived. 

PACS numbers: 42.50.Dv, 42.50.Lc, 42.50.Pq 



I. INTRODUCTION 

Methods of detection for optical fields such as the het- 
erodyne (and homodyne) [l[ enable measurements of the 
field amplitude (and quadratures) and are useful for de- 
tecting nonclassical light [5, 's'l . The heterodyne (or ho- 
modyne) detector is an example of general diffusive mea- 
surements d, Q , so-called because the statistical fluctu- 
ations in the current outputed by the measurement is 
driven by a Gaussian-white-noise process, a prime exam- 
ple of a diffusion process. This is the process responsi- 
ble for the broadening of a probability distribution de- 
scribed by the diffusive term appearing in Fokker-Planck 
equations @, 0]- The heterodyne technique, originally 
conceived for radio technology more than 100 years ago 
[1, Q, of which the homodyne is a variant, has now be- 
come an indispensable part of many quantum informa- 
tion processing applications. Some of these applications 
include quantum-state tomography Il0l- |l2l| , teleportation 
[T3l - [l6| . and state preparation 17[-|22|. It also enables 
tests of quantum mechanics [23l - 2a | . 

In this paper we are concerned with continuous mea- 
surements (also referred to as monitorings) of an open 
system in a vacuum bosonic environment (Fig.[T]). Specif- 
ically we are interested in how one can specify the class of 
all possible diffusive quantum measurements given how 
many inputs to the measuring device there are. There 
is already one solution to this problem [2^ [23|, a pa- 
rameterization that we shall call the U-rep. The U-rep 
parameterizes an arbitrary diffusive measurement by a 
square, real matrix U, with three constraints imposed on 
the blocks of U. We will present two more solutions with 
certain advantages over the U-rep: The M-rep — a non- 
square complex matrix M ; and the B-rep — a square ma- 



trix and two vectors. We will be comparing the different 
parameterizations in detail later but let us first say why 
one might consider these alternative parameterizations. 

A good reason, and also the more theoretical reason 
for considering the M-rep, is that aesthetically it is more 
attractive than the U-rep. As we will show, the M-rep 
has the ability to unify all diffusive measurements by a 
single equation whereas three would be required to define 
U [75| . The aesthetic advantange of the M-rep continues 
when we come to the quantum theory of multiple-input 
multiple-output Markovian feedback control based on dif- 
fusive measurements [28j . There, the theory is simpler 
when the equations are expressed in terms of the M-rep 
rather than the U-rep. Furthermore, the defining equa- 
tion for M has an intuitive interpretation which makes it 
easy to remember. 

Either the U-rep or the M-rep can be used when we 
know what the measurement is and would like to param- 
eterize it in order to model the state of the open system 
conditioned on the results of the monitoring. If on the 
other hand we were given either a U, or an M, and were 
asked to describe the measurement with actual optical 
elements then this is a much more difficult task for ar- 
bitrary measurements. Such a setting where one may 
be given an M without knowing its physical implemen- 
tation can in fact arise naturally in quantum feedback 
control [23, : As measurement is an inherent part of 
the feedback loop, a design of the feedback loop can thus 
incorporate a design of measurement and it is natural 
to ask what sort of measurement one should do in or- 
der to achieve a particular objective for the control. In 
the case of optimal control the objective would be to 
minimize a measurement-dependent cost function. The 
result of this optimization would be a matrix M*, which 
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one would then need to realize in the laboratory. That 
is, the theorist who now has obtained M* would like to 
inform his/her experimentalist colleague about how to 
construct the measurement. 

The above considerations motivate us to propose yet 
another parameterization of diffusive measurements; one 
which we call the block-rep, or B-rep for short. Un- 
like the M- and U-rep, the B-rep is a realization of the 
measurement in terms of beam-splitters, phase shifters 
and homodyne detectors. The B-rep is so-called because 
it parameterizes the diffusive measurement in terms of 
three distinct blocks, with each block corresponding to a 
set of parameters (beam-splitter transmission coefficients 
and phase shifts). Note that being given an M is equiv- 
alent to being given a POVM (positive-operator- valued- 
measure [7^). From this point of view the relationship 
between the M-rep and B-rep is thus one of a POVM and 
its realization. As we shall see, proving that an arbitrary 
POVM has a realization defined by the B-rep is much 
more difficult than translating a given B-rep into an M, 
and hence a POVM. Here we point out that this line of 
thought had in fact been applied to A^-port homodyne 
detection — a generalization of the standard homodyne 
measurement but a subset of all diffusive measurements, 
as early as 1987 3l|. A similar question was raised again 
in 1994, but for unitary operators instead of a POVM 
[33j . Our construction of the B-rep in fact relies on this 
resuh of Ref. [H. 

Our paper is organized as follows. In Sec.|II]the uncon- 
ditional dynamics of open quantum systems is reviewed 
and the concept of unravellings and its relation to con- 
tinuous measurements are also briefly sketched. As we 
frequently make use of vector-operators, henceforth re- 
ferred to as vops, a comprehensive review of their def- 
initions and algebra is provided in Appendix |Al Note 
that not all of the results in Appendix |X] will be used in 
this paper so we will refer to the ones that do appear. 
However, the reader who is also interested in the the- 
ory of multiple-input multiple-output quantum feedback 
control (for which the M-rep is applied to) [28] may find 
it worthwhile to embark on a fuller exposition of vop al- 
gebra. We also mention that for convenience we will not 
necessarily reflect the multi-component nature of vectors 
or vops in our language when they are referred to, such 
as in "the field a", or, "the current y", as opposed to 
using plurals as in "the fields a" or "the currents y". 
Our first results section begins with Sec. IHIl The key re- 
sults here are the definition of the M-rep and its relation 
with the U-rep. It should be noted that Sec. IHIl begins 
by postulating the most general diffusive stochastic mas- 
ter equation. The purpose of Sec. [Tl] is thus to help the 
reader gain the intuition required to make this leap. In 
Sec.|lV]we define the B-rep and discuss its relation with 
the M-rep. Here we conjecture a matrix decomposition 
of M in terms of the matrices in the B-rep triple. We 
then conclude in Sec. IVl 



II. DIFFUSIVE CONTINUOUS 
MEASUREMENTS IN THE SCHRODINGER 
PICTURE 

A. Unconditional System Dynamics 

A genera l op en quantum system with Markovian dy- 
namics [33, l35j | can be described by a master equation in 
the Lindblad form [H, 113 



hp - 



[H,p]+V[c\p=Cp 



(1) 



where H is Hermitian and c = (ci, C2, . . . , c^)^ (which 
we will refer to as a vop) is a vector of Lindblad opera- 
tors, all time-independent. See Appendix |X] for a com- 
prehensive review of definitions and conventions that we 
are adopting for vops. We have also defined 



2)[c] = ^I?[cfc], V[c\p = cpc^ 



\ c^cp - 



ipc+c. (2) 



k=l 



Note that instead of working in natural units by setting 
h — 1 we will work in units such that a factor of H ap- 
pears on the LHS of ([1]). This was done in Refs. 0, 
so we will keep their convention for ease of comparison 
with the results therein. The reason for redefining units 
this way is to keep track of the correspondence between 
quantum operators and their classical counterpart (the 
dynamical variables which they represent) and aspects 
of quantum control that are not present in the classical 
theory. Equation ([T]) can be derived by a unitary op- 
erator U{t,to) ac ting on the joint Hilbert space of the 
system and bath [33]. This describes the interaction be- 
tween the system and its environment, which we assume 
to be given by the Ito stochastic differential equation 



hdU{t,to) = 



'-iH -hc^c 



b^c - c%)U{t,to)dt , 



where dU{t, to) = U{t + dt, to) - U{t, to) and 
= (cj, C2, . . . , c\) ■ 



(3) 



(4) 



Here b represents a bosonic bath field shown in Fig.[T]and 
assumed to be in the vacuum state. External driving 
of the system is included in H. In this section we are 
working in the Schrodinger picture so operators are time- 
independent, equal to their initial value. 

It is well known that the evolution described by ([T]) 
can be decomposed into stochastic paths in Hilbert space, 
called quantum trajectories (a term coined by Carmichael 
[ssl who also referred to this decomposition as unravel- 
ling the master equation) by considering the evolution 
of p conditioned on the results of monitoring [H, |39| . 
We will also refer to the results of monitoring as a cur- 
rent, denoted by a vector y(i). For Markovian evolution 
it is sufficient to consider unravellings (the set of possi- 
ble quantum trajectories the system state may take over 
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FIG. 1: Schematic of some photoemissive system which cou- 
ples to the environment via some set of operators c. The time 
evolution of the system based on knowledge of the noisy sig- 
nal y also follows a noisy path, which can be described by a 
stochastic master equation. The stochastic master equation 
updates the system state continuously in time as the mea- 
surement record grows. If one ignores the measured result 
y then the system state evolves according to the determinis- 
tic master equation. This corresponds to averaging over the 
stochastic evolution. 



time) generated by measurements with either or both of 
two classes of noise 

0, 113: 1) A point (or jump) pro- 
cess also called a Poissonian random variable, or; 
2) a Wiener process also called Gaussian white noise. 
We will refer to this second case as a diffusive measure- 
ment. Physicists usually neglect the technical difference 
between an "unravelling" and a "measurement," and one 
often finds the two terms used interchangeably. We will 
do so as well in this paper for convenience. 



B. Conditional System Dynamics — Homodyne 
and Heterodyne Unravellings 

The most familiar forms of diffusive measurements in 
quantum optics are homodyne and heterodyne detection 
schemes of a single output field of some photoemissive 
source. Apart from being characterized by Gaussian 
white noise (as opposed to point-process noise) these 
measurements are also determined by physical properties 
such as the specific arrangement of linear optical elements 
and the specific parameter values at which they operate. 
Here we will first show how the master equation ([Ij is 
unravelled by the homodyne and heterodyne detection 
schemes. This will familiarize the reader with the nota- 
tion used for describing conditional evolution and build 
some intuition about how the measurement parameters 
enter into a quantum trajectory equation. For conve- 
nience we will set h = 1 for the rest of this section. 



1. The homodyne stochastic master equation 

We begin with the example of a single mode of an in- 
tracavity field. For simplicity we will assume the optical 
cavity to be one-ended with transmission coefficient 7 
and that the bath is in the vacuum state. In this case 
the master equation ([1]) has only one Lindblad term (i.e. 
L = 1), given by c = ^/ja, where a is the cavity mode 
annihilation operator. A homodyne measurement [l| of 
the field at the leaky port would then give us information 
about a particular quadrature of d. Such a continuous 
measurement would yield a result, for some state p{t), 
given by the Ito formula [1, (assuming a measure- 
ment efficiency of 77) 

y{t)dt^ yffj {c + c^){t)dt + dw{t) . (5) 

We have defined, for an arbitrary operator A, {A){t) = 
TT[Ap{t)], reserving angle brackets for quantum (opera- 
tor) averages. The term dw{t) is a Wiener increment. It 
is a Gaussian random variable defined by the mean 

E[dw{t)] = (6) 

and variance 

[dw(t)]^ = dt , (7) 

dw{t)dwit') = yt^t'. (8) 

We shall always denote an average taken with respect to 
a probability distribution by an E as in (|6]). 

From the measurement result y{t) the observer can 
update her prior state from p(t) to a posterior state 
Pyt + d,^) according to the stochastic master equation 

dpy, {t) = C p{t) dt + y/T] n [c] p{t) dw , (9) 

where C is given by ^ and the change in the state is 

dpyM = PyAt + dt)^ p{t) . (10) 

We have also defined (see also (|A65p in Appendix El) 

= is + sit -Tr (is + sit) S , (11) 

for any A and B. Here py^ {t + dt) is the state given a par- 
ticular realization of y{t). We have also written the time- 
dependence of y as a subscript. For clarity it is best to 
have this flexibility in writing the time-dependence. Av- 
eraging ^ over y{t) according to its actual distribution 
piijt) (a Gaussian with mean y/rj (^c + c'^'j and variance 
1/dt) simply returns the master equation ([T]). Note that 
a breve over a random variable denotes its realization. 

Equation ([9]) is a nonlinear stochastic differential equa- 
tion in the state. This would not have been the case if we 
did not insist on keeping the state normalized at all times 
and having the measurement results distributed accord- 
ing to its true statistics described by p{yt)- An alterna- 
tive theory is to assign an ostensible distribution postijjt) 
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to y{t) and just disregard the norm of the state. When 
Postivt) is a Gaussian with zero mean and variance 1/dt 
the state evolves according to a linear stochastic master 
equation given by 

dpy,{t) ^ Cp{t)dt + y{t)dtn[c\p{t) , (12) 

where the updated state is now unnormahzed 

dpyM=PyAt + dt)-p{t), (13) 

and we have defined a "hnear" version of (|lip 



n[A]B = AB + BA'' 



(14) 



The solutions to (112]) are called hnear quantum trajecto- 
ries (4^ . For hnear quantum trajectories to be equivalent 
to the standard theory (the case when the measurement 
results are distributed according to their true distribution 
piijt) and the state fohows a nonlinear stochastic master 
equation) the ostensible distribution must be such that 



p(yt) dyt = Tr[py^{t + dt)] postivt) dyt . (15) 

The form of the linear stochastic master equation is spe- 
cific to the form of post(i/t), so choosing a different form 
for post(yt) will result in a different linear stochastic mas- 
ter equation 

In the above we have assumed the prior state to be 
unconditioned (implied by (jlOp ). no matter how it is ob- 
tained. In general this need not be the case, and certainly 
will not be once the prior state is propagated for any finite 
time interval under (or ^ ) . Thus for generality and 
clarity we will introduce a Roman subscript c (for "condi- 
tioned") and write the linear stochastic master equation 
as 



where 



dpc{t) ^Cp{t) dt + y{t)dtn[c]p{t) 



dpc{t) = p^{t + dt) - p{t) 



(16) 



(17) 



Tacked onto Eqs. p6)) and (|T7l) is the system of nota- 
tion that p{t) denotes a state of knowledge that is ar- 
bitrarily conditioned on the measurement record during 
[0,i). Once p{t) is determined the precise conditioning 
of pc{t -\- dt) can then be expressed. We will inform the 
reader about the dependence of p{t) on the measurement 
record when it is necessary to assume a particular condi- 
tioning for p{t). 



This latter detection scheme of using two homodyne de- 
tectors is also referred to as a dual-homodyne measure- 
ment. It is a convenient way of understanding the het- 
erodyne measurement and is how we will think about 
it. That is, while we use the term "heterodyne" often, 
we will always realize it using a dual-homodyne scheme. 
This decomposition of a heterodyne measurement into 
two homodyne ones will also be of use later in Sec. IIVI 

If we now keep the above example of a leaky cavity and 
perform a heterodyne (or dual-homodyne) measurement 
of the output of the cavity then we obtain two currents, 
yi and 2/2, which we write as the components of a 2 x 1 
vector y. 



ydt = 



(c + ct) 



dt + dw . 



(18) 



where dw = {dwi,dw2)^ and dwi and dw2 are indepen- 
dent Wiener increments. 



dwi{t)dw2it') = yt,t'. 



(19) 



Conditioning the system state on (|T8)) then leads to the 
heterodyne stochastic master equation 



dpc{t) =Cp{t)dt 

^ 'H[c\p{t)dwi + 'H[~ic\p{t)dw2 ■ 

(20) 

We can write this more compactly by defining the row 
vector 



M = a/^ (1,^) 



(21) 



Equation (f20|) is then equivalent to 

dpc(0 ='Cp(i)dt + dw^H[Mtc]p(i) . (22) 



Note here the use of the vop- valued superoperator (|A66p . 
defined for any A and B as 

n[k\B = AB + BA^ -TT:[kB + BA^]B , (23) 

where = (A^)^. The heterodyne measurement illus- 
trates how y can be a two-component vector even though 
L = 1 so there is only one output field feout- 



III. THE M-REPRESENTATION 



2. The heterodyne stochastic master equation 

A heterodyne measurement of efficiency 77 can be 
shown to be formally equivalent to two simultaneous 
measurements of orthogonal quadratures of the signal 
field by two homodyne detectors each with efficiency 77/2. 



In this section we show how a single complex matrix 
can capture any such properties of diffusive measure- 
ments of an arbitrary number of system outputs in the 
Schrodinger picture. In Sec. IIVI a second and physically 
intuitive method of parameterizing diffusive unravellings 
is formulated in the Heisenberg picture. 
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A. Measurement Statistics 

Based on the foregoing examples of homodyne and het- 
erodyne unravehings in Sec. Ill Bi we propose the following 
form of a general diffusive stochastic master equation 

hdpcit) ^ Cp{t)dt + dvf^'H[M'^c]p{t) , (24) 

conditioned on the 2L x 1 real vector 



y{t) dt^E [y{t)] dt + dw{t) , 



(25) 



which represents the measured current. The length of 
y is motivated by the example of heterodyne detection 
in Sec. Ill B 21 and the fact that here we are allowing L 
components in bout- The measurement is a noisy pro- 
cess. The noise in the measurement is given by a vector 
of independent Wiener increments dw (also 2L x 1 by 
default). The vector d'w{t) therefore has zero mean 



E[dw{t)] =0 Vi 



(26) 



and correlation matrices 



dw{t)dw'^{t) = l^^dt , 
dw(i)dw^(t') = yt^t' 



(27) 
(28) 

In ([24]) we have used ([23]) and defined M to be an L x 2L 
complex matrix that parameterizes the unravelling (77| . 
The constraints which define M will be determined in 
the following sections, here we will obtain the form of 
the measured current consistent with ((24)) . 

The precise of form of E[y{t)] can in fact be derived 
from the theory of linear quantum trajectories which al- 
lows us to calculate the actual statistics of the measured 
current by using an ostensible distribution in the follow- 
ing way: Here we are only concerned with the average of 
yt and this is 

E[y{t)]dt = E[y(t)Tr{p(i) + dp^)}] dt , (29) 

where p{t) is normalized and dpc(t) is given by 

hdpcit) ^ Cp{t)dt + y'^{t)n[M''c]p{t)dt. (30) 

Here we have used (fT5|) to rewrite the average on the LHS 
as an average with respect to the ostensible distribution 
(denoted by an overbar on E) 



post {ft 



dt 



^ exp[-^yjytdt] . (31) 



This is a Gaussian with mean zero and covariance 

E[y{t)y^(t)]=hjdt. (32) 

Defining a linear version of (|23)) , 

H[A]B = AB + BA^ , (33) 



we then obtain, on substituting ((30)) into ((29)) . 

E[y{t)]dt = h-^E[y{t)dt y^{t)dtTT{H[M'< c\p{t)}] 

(34) 

= h-^TT{M^cp{t) + p{t)M^c^} 



dt 



= — (Mtc + M^c*) 



(35) 



We have used the facts that Cp is traceless and that 
E[y] = in ([M]). Later (Sec. |lV]) we wih be consider- 
ing a different representation of diffusive measurements. 
In order to distinguish between different representations 
we will use a subscript on the current. Substituting p5p 
back into (1^51) we will thus write 



hyu dt = {M^c 



re* 



)dt + hdw 



(36) 



It is clear from p5l) that M determines what property 
of the system gets measured, represented by some Her- 
mitian operator /(c), and also the measurment statistics 
(i.e. the statistics of yt) for a measurement of /(c). That 
M defines the diffusive measurement could have been an- 
ticipated from ((M)) since M appears in the H superop- 
erator which stems from considering measurement. Its 
appearance in ((^^ also means that not every M will be 
a valid parameterization of a diffusive measurement. In 
Sec. IIIIBI we find the necessary and sufficient condition 
for M to be valid. 

The careful reader who is familiar with the quantum 
theory of indirect measurements [4, 43] will know that the 
conditioning for the state is in general just a label which 
allows us to distinguish between measurement outcomes, 
and as such it does not have to be real. Our choice of 
conditioning the state on a real vector is motivated by its 
use in control theory in which the prevalent treatments 
regard y as a real process fil ]. 

The assumption that p{t) is unnormalized can be re- 
laxed. Say the measurement began at time and that 
Pc{t) is obtained by evolving p(0) via ((30)) to time t. In 
this case the correct modification to the above calcula- 
tion is to divide the trace in ((^ by Tr{pc(i)}- Keeping 
the same ostensible distribution as ((?T)) we get 



E[yt|y 



[0,t) 



Idt^E 



yt 



TT{dp,{t)} 



Tr{pc(t)} 



dt 



(37) 



where on the RHS the state pc is conditioned on the 
measurement record 



y[o.o = {y(s) I < s < t} . 



(38) 



The end result from calculating ((57)) is to replace the 
unconditioned averages in ((35p by conditional ones, 

(39) 



(c)^(c)e = Tr[cpe(i)] 
as one would have guessed, where 

Pcit) 



Pc{t) = 



Tr[Pc(i)] 



(40) 
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B. The Set of Allowed M 

1. Necessary condition for M 

If M is to be a valid parameterization of a diffusive 
quantuni measurement it must be such that (I24p evolves 
a valid state at time t to another valid state p{t + dt) 
for all t. Here we remind the reader that a valid state is 
represented by an Hermitian operator that is normalized 
and positive. From this it follows that its eigenvalues 
all lie in the interval [0,1], which we can write as an 
operator- inequality < p < 1. The positivity condition 
is the only nontrivial criterion because Cp and T-Lp are 
each Hermitian and traceless so (fM)) will always preserve 
Hermiticity and normalization for any M. Since these 
conditions must hold for any measurement process, they 
are necessary for a given M to be valid, and we establish 
the following implication 

M is valid ^ M : < p{t) + dpc{t) < 1 . (41) 

where p(t) is understood to be any valid state and dpc{t) 
is given by It is simple (given dpc{t) is traceless) to 
see that 

< p{t) + dp,{t) < i =^ Tr{[p(t) + dp,{t)f} < 1 . 

(42) 

Without loss of generality we can assume that p{t) is a 
pure state, because of the convexity of Tr[p^] in p. We 
will show 

Tr{[p{t) + dpc(t)]'} < 1 ^ MMVri e io , (43) 
where Sj is 

jO = {H-diag(»7)|Vfc, 77fee [0,1]} , (44) 

thereby establishing the RHS of the equivalence as 
a necessary condition for M to be valid. 

The equivalence ()43p can be proven by first proving 
the following lemma; 

Tr{[p(i)+dpe(t)]'} = 1 ^ MMt/r^-I, , (45) 

where 1^ is the L x L identity matrix. It can then be 
shown (see Appendix [B|) that 

Tr{[p(i) + dpc(t)]'} = 1 + ^ tr[ff (MMt/^ - I^)^] dt , 

(46) 

where H is the Hermitian, positive-semidefinite matrix 
(not to be confused with the Hamiltonian H, which is an 
operator) 

H={{c-{c))\c-{c)f) . (47) 

Since p{t) is an arbitrary pure state, H can be assumed 
to be strictly positive (i.e positive definite). Thus pS)) 



will be 1 if and only if MM^ /h = 1^ . This completes the 
proof of (PS]) . a necessary condition on M in the case of 
efficient monitoring (i.e. monitoring which preserves the 
purity of the state). 

To prove (|43|) we must consider the case when the pu- 
rity at time t + dt drops to below 1. Given that we had 
a pure state at time t this is possible if and only if we 
had inefficient monitoring. For a single decay channel we 
can bundle the sum of all losses [tI] into a single parame- 
ter 77, and consider only a fraction rj (between and 1) of 
the system operator c to be measured perfectly. This can 
by modelled by introducing an imaginary beam-splitter 
with transmission coefficient rj in the path of the decay 
channel followed by an ideal detector [4^. When mul- 
tiple decay channels are present we simply repeat this 
setup for each channel. This motivates us to rewrite the 
master equation ([T]) as 

hp^^i[H,p]+V[y/T~Hc]p + V[VHc]p , (48) 

where H e [recall (|44l) ]. and unravel the last term with 
unit detection efficiency. Such an unravelling is defined 
by the stochastic master equation 

hdpc = Cpdt + dw^n[M'''c]p , (49) 

where c' = \/H c and (H^ is conditioned on the current 

y^^, dt=j (M'^c' + M'^c'*) + dw . (50) 

Since (|49)) and (|50l) describe efficient monitoring, from 
the lemma of (|45|) . M' must satisfy 

M'M'V^=Ii,- (51) 

The unravelling defined by (|49|) . ([50]) . and (|5T|) is equiv- 
alent to dMl), (ESI), and (I2S1), but the latter make the 
measurement efficiency H explicit. The set of quantum 
trajectories generated by (|24|) must therefore be the same 
as the set generated by (|49p . This will be the case if and 
only if the two equations have the same stochastic term, 
i.e. 

dw^-H[Mtc]p = dw^-H[M'^cV. (52) 

This is true for any p and c if and only if M = \/H M'. 
This gives 

MM^ = H. (53) 

We have arrived at (|53|) by considering inefficient detec- 
tion. However there are other properties of the measure- 
ment that one would like to capture with M so it would 
seem that ([53]) is necessary but not sufficient. As we show 
next, ()53p is in fact, surprisingly, a sufficient constraint 
on M for it to be a valid parameterization of general dif- 
fusive measurements. 
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2. Sufficient condition for M 



To show that ((53|) is a sufficient condition we wih con- 
nect M to another parameterization U, which is a differ- 
ent way of representing diffusive measurements and has 
conditions that have previously been shown to be neces- 
sary and sufficient [2^, [13] ■ In this paper we introduce 
an elegant way to connect M and U by considering a 
generalized diffusion operator which we denote by D. If 
the system state is in an A^-dimensional Hilbert space H, 
then D is an operator in H (g) H, defined by 



D dt — dpc <8) dpc 



(54) 



This is the mathematical object which characterizes the 
set of all equivalent representations (i.e. M-matrices) of 
a given unravelling as we now explain. 

Stochastic paths of the quantum state itself are rather 
abstract but we can make the trajectories more concrete 
by considering 



Pc = Tr[epc] , 



(55) 



where e = (ei, 62, ... , ejv^)^ is an operator-basis for all 
linear Hermitian operators. Thus pc is a stochastic pro- 
cess in , satisfying 

dpc = Adt + Bdvv . (56) 

Note that A = Tr[e£/9] is a vector while 

B = Tr{eH[c^M]p} (57) 

is a matrix. Recall that £ and H[c^M] are defined by 
([!]) and (|23|) respectively. The diffusion matrix for pc is 
given by 



Ddt^ dpc dpj ^ BB^dt 



(58) 



We can see how D arises by rewriting the RHS of the first 
equality in ([58)) . 



Ddt = Ti-[edpc] Tr[e^(ipc] 

= Tr[(edpc)«)(e^dpc)] = Tr[(e (g) e^) 6] . (59) 

Thus D is an operator- valued diffusion coefficient whose 
trace against e®e^ gives us the diffusion of the more tan- 
gible process ((55|) . Turning pc, a stochastic process in H, 
into Pc an equivalent stochastic process in M.^ , allows 
us to grasp the abstract diffusion operator D by using 
well known properties of classical stochastic processes; 
namely that the diffusion of Pc is characterized by the 
matrix Z?, and that for a given D there are many ma- 
trices B such that BB^ = D. The equivalence between 
Pc and Pc and more specifically between B and M, and 
also between D and D (as given by ((57)) and (l59t ). mean 
that D characterizes the diffusion of pc just as D charac- 
terizes the diffusion of Pc- There will be many matrices 



M that generate the same D (through the nonlinear term 
H[c^M]p) just as there are more than one B correspond- 
ing to a given D. 

We mentioned earlier in Sec. Ill Al that an unravelling 
can be defined as the set of solutions to the stochastic 
master equation and this set will correspond to some M. 
In general the set of solutions will change everytime M is 
changed but it may also be possible to have two different 
choices of M, say Mi and M2 which generate the same 
solution set. In this case we say that Mi and M2 are 
equivalent representations of the same unravelling. Let 
us substitute ([M)) into (IM)) and define 



T = 



Ti 
T2 



_ /3?[M] 



3?[M] 



(60) 



where Ti and T2 are each real, and L x 2L so that T 
is a real 2L x 2L matrix. We have defined the real and 
imaginary parts of an arbitrary complex matrix A by 

^[A] = ^ {A + A*) , ^[A] = ^ {A - A*) . (61) 

We can then write M as 

M = Ti+iT2. (62) 
The vop M^c can be rewritten as 

M^c = T^c& , (63) 

where = ( c^, —ic^). Using this, ([M)) is 

Ddt ==-H[cjTdw]p®-H[dw^T^c&]p 
= H[c1]pTdw^ dw^T^H[c&]p 
^hH[cl]p (g)UH[c&,]pdt . (64) 

In the last line we have defined a 2L x 2i real matrix U 
by 



TT' =nx} = h 



Uii U12 
U21 U22 



(65) 



where Umn are L x L blocks of U. Equation (|64l) tells 
us that given the set of Lindblad operators c, U is what 
uniquely specifies D, and therefore the unravelling. That 
U defines the unravelling was independently formulated 
in Rcf. [23 for the case of pure-state trajectories, and 
later generalized to include non-unit detection efficiency 
in Ref. [27|. The necessary and sufficient conditions for 
a U to be valid (in the same sense as (|1T|) ) are 

U > , (66) 
Uii + U22 e -So , (67) 
U12 = U21 . (68) 

We will denote the set of valid U-matrices as 
Uii U12 



U21 U2 



U > 0, Uii + U22 e .$3, U12 = U21 



(69) 
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The question is now whether any U derived from ((53)) 
is valid. Since T is a real matrix, the definition (j65p 
ensures (|66l) for any T. Using (|65l) conditions (f67| and 
become, respectively. 



which is consistent with (l69l). 



TiT^+TsTjeiO, 
TiTj = T2T7 . 



(70) 
(71) 



This defines the set of valid T-matrices 



Ti 

T2 



(72) 



It is easy to see, on substituting (|62|) into (l53l) . that ([70)1 
and ([7T|) are satisfied and thereby showing the sufficiency 
(and in fact necessity) of (f53| for M to be a valid repre- 
sentation of an arbitrary diffusive unravelling: 



U(M) eU ^ Mem 
where DJl is the set of valid M-matrices 
m={M\ MMt e io} . 



(73) 



(74) 



Note that even though we have imported a previous re- 
sult which proves altogether the sufficiency and necessity 
of (|53p . our analysis in Sec. IIII B ll remains instrumental; 
since without it there is no reason to consider the set DJl. 
For ease of reference we call the members of DJl "M-reps" , 
and similarly the members of ii "U-reps" . 



C. Comparison between and il 

Here we discuss some general features of M-reps and 
compare them with U-reps. To do so we first recall the 
essential features of the U-reps from Re fs. p6l . [27j . There, 
Wiseman and co-workers derived a stochastic master 
equation of the form 



dpc = Cp{t)dt + n[d-i^c]p{t) 



(75) 



(with C given by ([T])) that conditioned on an L x 1 com- 
plex current 



3dt = {Uc + Yc}) dt + dz , 



(76) 



where dz (also L x 1 by default) is a complex white-noise 
increment defined by E[(iz] = and the correlations 



dz d-z.^ = hUdt , dz dz^ = HY dt 



(77) 



The matrix H is a member of ^, while Y is complex 
and only has to be symmetric. The matrix U, or the 
"unravelling matrix" as it was originally called, is defined 
as the correlations of (jR[dz^],'^[dz^f) /^/h . This can be 
shown to be 



U = 



1 



H + 5R[Y] 



5[Y] 

H - my] 



(78) 



The first noteworthy feature of M-reps is that they cap- 
ture physically valid measurements by a single equation, 
(l53l) . whereas U-reps require three, namely (|66)) - (l68t . 
This point makes our theory of M-reps rather elegant. 
If however we compare the number of real parameters 
required to specify M to that of U, then we find that 
U-reps are a more efficient parameterization: We know 
directly from ()78p that the number of real parameters in 
U is + 2L, being the total number of real parame- 
ters required to specify H and Y. To find the number of 
real parameters required to specify M, we note that M 
is L X 2i, so the total number of real parameters in M 
is 4L^ . But ([5^ means that — L of these parameters 
are redundant (since it imposes I? — L independent real 
conditions on the elements of M). Subtracting this re- 
dundancy from leaves 3L^ -I- L. Thus M-reps require 
21? — L more parameters than U-reps. 

The second attractive feature of M-reps is that M can 
be understood as the generalization of -v/R , encapsulat- 
ing more information than just measurement efficiency. 
This was in fact part of our motivation for developing 
M-reps and postulating ([M]). as stated at the start of 
Sec. IIII Al So in hindsight, the fact that M turns out to 
be the matrix square root of H, i.e. constrained by (j53p . 
seems "natural" given that we were thinking of gener- 
alizing to -v/H , and then from \/H to M from the 
start. The interpretation of U-reps on the other hand 
was derived (summarized in ([75)) - (|78)) ) as the correlation 
between the real and imaginary parts of a complex mea- 
surement output. This has the advantage that it reflects 
directly the quantum trajectory diffusion as was evident 
in (|64|. Equations (|66|) - (l68l) are simply constraints on 
the allowed correlations of dz. 

Finally, we mentioned before (two paragraphs under 
P6p) that we prefer to work with a real current, namely 
Ym, for the purpose of feedback control even though 
quantum measurement theory does not impose such a 
restriction. However, by allowing the system state to be 
conditioned on a real vector, yM now has the capacity to 
be an observable (i.e. correspond to an actual physical 
process) rather than a mere dummy variable for distin- 
guishing the different measurement outcomes. Of course, 
one may also prefer a real current over a complex one sim- 
ply because it is real. Motivated by exactly these reasons 
Ref. 2Ji defines a real current from J as 



(79) 



where T"*" is the Moore-Penrose inverse of T [i^. Sub- 
stituting ((76)) into (|79| and using the properties of T+ 
it can be shown that (l79)) simplifies to the form of (|36)) 
(recall (|62))). Thus the U-rep achieves the same current 
as the M-rep but less directly by deriving y via J. 
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D. Relations between T, 11, and 3K 

The relation between 9Jl and 2 is straightforward, 
given by ([50)1 and (p^ . From these it should be clear 
that DJt maps one-to-one and onto 1. We said above 
that M-reps are parameterized by + L real numbers, 
therefore T also has the same number of real parameters. 
This number is 2L^ — L more than the number of real 
parameters in the U-rep. We can see where the extra pa- 
rameters lie by returning to (j65|) which defines T as the 
matrix square root of U. We can always write, using the 
polar decompostion, T = -\/TT^/?i = -s/U where VU 
is the positive matrix square root of U [79j and is some 
2L X 2L orthogonal matrix which may or may not be 
unique given T (see following paragraph). The number 
of parameters in \/U is + 2L, inherited from U, while 
the number of parameters in is 2L^ — L [sO] . Adding 
the number of parameters in -s/U and gives 3L^ -1- L. 
The extra parameters thus lie in 0, leading us to define 
the set of pairs (not to be confused with a block matrix) 

iio = {(%/U,0)|Ueil,0^ = 0-i} , (80) 

where is understood to have dimensions 2L x 2L. 

We can divide the set 1 into two subsets: a subset 
of invertible Ts and a subset of non-invertible Ts. The 
subset of Ts that are invertible maps one-to-one and 
onto ilo, since in this case the polar decomposition of 
T is uniqu e, i.e. given any T, both entries of the pair 
(■y/TT^/?i , 0) e ilo are uniquely determined by T. The 
subset of Ts that are not invertible maps infinitely-many- 
to-one and onto a subset of ilo- In this case the first en- 
try of the pair {y/TT^ /h , 0) will still be uniquely deter- 
mined by T but corresponding to this T there are an in- 
finite number of orthogonal matrices 0. We depict these 
relations in Fig. [2] along with another parameterization 
*B to be introduced in Sec. IIVI Note the cardinality of 
the non-invertible subset of T is less than the cardinality 
of T. That is, the former is a subset of measure zero. 

E. The M-representation in the Heisenberg Picture 

We now consider diffusive quantum measurements in 
the Heisenberg picture where operators evolve and states 
are static, equal to their initial value. Moving from the 
Schrodinger picture to the Heisenberg picture requires 
two conceptual changes which we summarize below. In 
Sec. IIVI we will formulate another representation of diffu- 
sive measurements in the Heisenberg picture so the dis- 
cussion here apply to Sec.|lV]as well. 

The bath field bin is now a vop- valued white-noise pro- 
cess. As such, dBin(t) = hin(t)dt satisfies Ito rules anal- 
ogous to (I2ni)-(I151) g3,l3. That is 

((iBi„(i)) - , (81) 
where the angle brackets denote a quantum average and 
dBinit)dBlit)^hhdt , (82) 




FIG. 2; Multivalued mappings between the different param- 
eterizations of diffusive quantum measurements. The lines 
have been drawn tangent to the circles to denote that one set 
is mapped onto the other, i.e. each member of one set has a 
corresponding member in the set to which they are connected 
to. Note also that we have defined a subset of T whose mem- 
bers are non-invertible. This is a set of measure zero and 
maps to subset of Uo (shown as an unshaded dot), also a set 
of measure zero. 



with all other second-order moments negligible. Note 
that 

[dBi„(t),dBf„(t')l =0 Vt^t' , (83) 

where we have used the mop-bracket (|A3ip . defined for 
any two vops A and B as 

[A, B] = AB^ - (BA^)^ . (84) 

Thus all other second-order moments with unequal times 
also vanish. As a consequence of the singularity of bin 
it is necessary to distinguish the bath field before and 
after its interaction with the system. This gives rise to 
an output field (Fig. [3]) 

dBout(t) ^uHt + dt, t) dBin{t) U{t + dt, t) (85) 

where U{t + dt, t) is defined by ^ but with b replaced 
by bin. From this one can show that 

dBout(t) = c{t) dt + dBin{t) . (86) 

This input-output relation and is what allows us to relate 
our measurement performed on the bath back to proper- 
ties of the system |48i, i51i] . 

We next turn the measured current y(<) into a vop 
[s^ - fs^ . Note that this does not mean we are taking 
the measured current to be a quantum object. The cur- 
rent is interpreted as a classical quantity but it is rep- 
resented by a vop in the theory. The classicality of the 
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FIG. 3: Quantum entanglement and measurement are two 
processes which have conventionally been described in the 
Schrodinger picture. When we move to the Heisenberg pic- 
ture we move to a difTerent mindset, where the system-bath 
entanglement and the measurement process are all captured 
by time-dependent operators, given by (|86|l [iol . [50| and (|87|l . 

measured current can be formally captured by the "self- 
nondemolition" property of the vop y(t) |55l - t57| . 

[y(i),y(t')l =0, Vt,t'. (87) 

It should also correspond to an observable process so y{t) 
must be Hermitian. In brief, quantum measurement the- 
ory in the Heisenberg picture just assigns an operator 
y{t) (a quantum stochastic process) which reproduces 
the correct measurement statistics (Le. the statistics of 
y{t), a classical stochastic process) Q. 

If a vop y is to represent the measured current y then 
it must be such that the statistics of y are reproduced, 
namely, 

(y(i))=E[y(t)] , (88) 

and 

[m dt] [m dt] ^ = [y{t) dt] [y{t) dt] ^1 . (89) 

Physically the measurement is never performed directly 
on the system but rather on the bath which involves de- 
tecting the output fields. Therefore y should depend on 
bout- It thus makes sense to write 

hyu dt ^ M^Bout + M'^ dBl^^ + h dvM ■ (90) 

We have defined the noise increment dvM to be 

dvM ^Idij + L* dlji , (91) 

where L is 2L x 2L, to be chosen so that ([55| and ([M)) 
are satisfied. The vop dXJ is a free vacuum field just like 
dBiii except that it is 2L x 1 so the only non-vanishing 
moment is 

dtS{t)dlj''{t) = hUr^dt . (92) 

Note that this means ((88)) is automatically satisfied. En- 
forcing we obtain 

/i^LL^ = fiU - M^M . (93) 



For convenience we will define Z = filai — M^M. We can 
thus choose fi. L to be the positive square root of the RHS 
of and write 

hl^VZ . (94) 

Note that dvM is not a Wiener increment (it has corre- 
lations given by Zdt/h) but we can write ym as 

hyu dt = (Mtc + M^c*)dt + h [rfv,n]M (95) 

where 

fi [dvmU = MtdBin + M^rfB*„ + hLdt + hL*dt^ 

(96) 

is a Wiener increment (when divided by h), i.e. 

MVm(t)]MMVm(0]M =l2Ldt , (97) 

[rfv^(i)]M [dv„,(i')]S = yt^t'. (98) 

We can show that yM satisfies the self-nondemolition 
property ([57)1 by first substituting either (^0]) . or ([^5)) 
into (|87)) and then considering the two cases t — t' and 
t t' separately. It will be quicker to use (|95l) and write 
(1571) in terms vacuum increments. Doing so gives 

[y{t)dt,y{t')dt] = [ [dv„(t)]M, Mv,„(t')]Ml = , (99) 

by virtue of ^ and ((55)) . 

Often one would measure, in the laboratory, the au- 
tocorrelation function of the photocurrent and its spec- 
trum (which is the Fourier transform of the autocorrela- 
tion function at steady-state), hence we now consider the 
autocorrelation function of yM- This is a two-time aver- 
age which can be calculated by using the vop quantum 
regression formula (jASSP . The derivation of the autocor- 
relation of yM follows the same method as in Ref. [l^ 
which derives the same quantity but with feedback. The 
only difference here is that these feedback terms would 
not appear so we will refer the reader to Ref. [1^ for de- 
tails. However we make the following heuristic argument 
which leads us to the correct form for the correlation 
function: First generalizing the correlation function of 
the homodyne current when L = 1 found previously in 
Ref. to the case of arbitrary L we obtain, 

(yM(t)y^(i + r)) 

= (^Tr|(VHc + \/Hc*) 

X e^'[c^V^ p{t)+p{t)c^VH]}Y + h^l,^6{T) , 

(100) 

where H € was defined in (|44| and C is given by ([ij. 
The action of the superoperator and the trace are defined 
in (jAlSp and (IA42P respectively. We said in Sec. IIIICI 
that M is a generalization of \/H . So the next logical step 
would be to replace -s/R by M. However this would not 
keep the Hermiticity of the terms standing on either side 
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of exp(£r), so the correct replacements in (jlOOp should 
be 



\/Hc — > M^c, VHc^ — > M^c* . I 
We then obtain 
^'(yM(t) YMit + r)) 

M*/9(<)+p(i)ctM]} 



101) 



|^Tr|(M^c + 
+ ^'U<5(r) 



(102) 



This is a generalization of the homodyne autocorrelation 
function found in Ref. [ssj ] to arbitrary diffusive measure- 
ments. 



IV. THE BLOCK-REPRESENTATION 

In this section we parameterize diffusive measurements 
by starting in the Heisenberg picture. The advantage of 
this is that it allows us to model the measurement by 
thinking directly about the physical transformations bout 
must undergo to produce a suitable current y. We show 
that the measurement block in Fig. |3] can be decomposed 
in a way shown in Fig. |4l 



A. Definition 

The blocks of Fig. |4]are as follows: 

1) Perhaps the first property of measurements that 
comes to one's mind is detection inefficiencies. For a 
single output this can be modelled by placing a beam- 
splitter with transmission coefficient rj in the path of 6out- 
The other input to the beam-splitter is not excited (i.e. 
no photons) so it is simply a field v in the vacuum state. 
We will choose a phase convention for the beam-splitter 
such that its output is 



y/rj bout + Vv V 



(103) 



where fj = 1 — rj. Note that we are retaining only one of 
the beam-splitter outputs since this is all that is neces- 
sary to model detector efficiencies. When multiple out- 
puts from the system are present (jl03p simply generalizes 
to 



H bo 



(104) 



where H = Ij, — H . This is the first block in Fig. 21 

2) Next, one is free to mix and phase shift the different 
components of bout to obtain an arbitrary linear combi- 
nation of them. Hence we multiply Sjn by an L x L matrix 
to get 



Sout 



Ssi, 



(105) 



We assume this mixing of the components of bout con- 
serves total boson number so that Sjj^Siu = Sout^out- This 
constrains S to be a unitary matrix. Note that S can be 
explicitly constructed out of L{L — l)/2 beam-splitters, 
L{L-\-l) /2 phase-shifters in a triangular arrangement js^l 
as shown schematically in Fig. |4l This scheme assumes 
that each component of bout have the same polarization, 
transverse mode structure, and mean frequency. If this 
were not the case then the appropriate corrctions would 
have to be made prior to entering the array in Fig. |4l 

3) We then split each component of Sout into quadra- 
tures. The process represented by K in Fig.|4]is shown in 
detail in Fig. [S] In general the intensities of the quadra- 
tures need not be equal so each beam-splitter in Fig. [5] 
has a different transmission coefficient which we denote 
by Ok (fc==l,2,...,L). This gives 



K 



Sout 



(106) 



If we_ define Q = diag(0) where 6k S [0, 1] for every k, 
and Q = li — Q, then K can be written explicitly as 



(107) 



Note that K is also a unitary matrix since dM 




t^out 



4) Individual homodyne detectors are then placed in 
the path of each component of d. The detectors are 
"tuned" to give a current that is normalized to the lo- 
cal oscillator amplitude that is inside it. Each homodyne 
detector is set to measure Re[c?fe]. This final current can 
therefore be defined to be 



y = 



d + d* 



(108) 



We see then, that once t], 9, and S are specified, it seems 
plausible that steps one to four above, by construction, 
will simulate an arbitrary diffusive measurement. Group- 
ing these parameters into a triple B = {rj, S, 6) we define 
the block-representation (or B-rcp for short) of diffusive 
measurements by 

{(»7,S,0)|Vfc, r]k,Ok e [0,1], St = S-i} . (109) 



5) We shall henceforth label the B-rep current in (jlOSp 
as yb . Once ye has been obtained one is free to do some 
post-processsing on this current. This will only be so if 
is an orthogonal matrix and the post-processing produces 
O^ye- We will have more to say about this in Sec. II V Cl 



B. Tiie Current Vector-operator 

Working out the input-output relation for each block 
in Fig. 21 we arrive at an equation for y that is a function 
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m 

4 



VL 



O y 



FIG. 4: Decomposition of the measurement box in Fig. [3] into blocks. Vacuum inputs are denoted by orange arrows. We are 
ignoring half of the output modes of H. Note that H and K do not mix the input modes (not counting vacuum inputs) i.e. each 
output mode of either H or K will not depend on more than one signal mode (mode coming from the left) and one vacuum 
mode (mode coming from the top). We represent this in the diagram by using highly asymmetric rectangulars. Each output 
of K then enters a homodyne detector, set to measure the input quadrature as defined in (|108p (also see Fig. [5} . 



(«out) 



Meaturemenb of 

VA-e {1.2... ..L} 




where 



Measurementtof 3?[r//,. 

VA-e {L + l.L + 2 2L 



FIG. 5: Expanded diagram of K. This prepares the compo- 
nents of Sout into quadratures, ready to be measured by units 
of homodyne detectors. 



of B = (t7, S, 6). Let us denote this dependence by using 
a subscript on y 



hyBdt = ^dBont + ^*dBl^^ + hdvB , (HO) 



J = Vh 



Q SVH 



(111) 



iy/QSVH J ' 

and hdvB is given by 

hdvB = ydV + y*dY^ +AdA + A*dAK (112) 

The noise increments c?V and dA are mutually uncorre- 
lated quantum Wiener increments. The coefficient ma- 
trices in (|112p are defined by 

V.^l-i^^^l, ,113) 

(114) 



Like the current in the M-rep, dvB represents correlated 
noise but we can write 





( ^qsTh 








( Vq ] 


7n 


{ -iVQ / ■ 



(115) 



hyB = J c + J* + [dvinle 
where 

[rfVm]B = JdBin + J*dBf„ 

+ VdV + V*dV* + AdA + A*dA* (116) 
has a mean of zero 

([dv„,(i)]B) = 0, (117) 
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and the only non-vanishing second-order moment 

[dv„,(t)]B [dvn,{t)]l = l^dt . (118) 

These can be understood physically by considering a vac- 
uum input, i.e. bout — bin- None of the transformations 
applied to the vacuum inputs in Fig. |4] alter their statis- 
tics, only add one vacuum field to another, which is again 
vacuum. We can also verify this intuition mathematically 
by using ((TTTt . pl4t . and (fTT3)) to derive (ITTSl) . Ah van- 
ishing statistics should be obvious from the properties 
of the vacuum inputs. One may also verify that (|87p is 
satisfied by ye- 

1. Example — dual-homodyne detection 

To illustrate how the B-rep works we will consider 
the familiar example of ideal heterodyne detection from 
quantum optics in the case of L = 1 , for which H = ry = 1 , 
S = e*"^, and Q = 9 = 1/2. Since the measurement has 
unit detection efficiency, there is no vacuum field v. It is 
simple to see that 

which transforms Sout = e"^6out and d into d = {di, ^2)^- 
From d the final output ye can be seen to be 

\y2j V2h\Xo,,t{(b + 7T/2)+^[a]J'^ ' 
where we have defined the the quadrature operator 

^out(<^) = {bout e'* + foLt e-*^)/2 . (121) 

C. Relations between S, 9Jt, and il 

Now that we have the M-re p, B -rep, and also the U-rep 
found previously in Refs. [26ll27l|. it would be instructive 
to study how one may convert from one representation 
to another. We can first derive a map *B — > 9Jt by not- 
ing that if the B-rep is to realize an arbitrary diffusive 
measurement then it must be the case that 

(ye) = (ym) . (122) 

Note that this is all we need since the second-order mo- 
ments have already been made equal. From (|90)) . (jllOl) . 
and (|llip we obtain 

M = jt = A/MTS^(VQ, -iv^O) • (123) 

One can verify from this that MM^ = hV\, as required 
if M is to be a valid representation of the measurement. 
Taking the real and imaginary parts of (|123p we get, 

V-VH3[ST]VQ -\/QS[ST]yH; 



From T we can construct U using the definition (|65l) . 
This gives 

Uii = VH 5R[S^]Q 5R[S] Vh + \/H 3[S^]Q 3[S] Vh , 

(125) 

U12 = - \/H 5R[S^]Q 5[S] %/H + \/H 5[ST]Q 5R[S] , 

(126) 

U21 = - \/H 3:[S^]Q 3f?[S] \/H + \/H 5R[S^]Q 3[S] \/H , 

(127) 

U22 = Vh 3[S^]Q 3[S] -I- yR 3?[S^]Q 5R[S] \/H . 

(128) 

By construction this U is positive-semidefinitc. By noting 
that if S is unitary, 

3fi[S^] 5R[S] + 3[S^] 5[S] = , (129) 
3fi[S^]3?[S] = 3?[S^]5R[S] , (130) 

we can verify that the remaining conditions for U to be a 
valid parameterization, ([57)) and are indeed satisfied 

by (in5|) -(fna. 

Equation (jl23|) says that given a B we can always find 
an M . But if we are given a valid M can we always find 
a B-rep? To answer this question we begin by counting 
the number of parameters required to specify B. The total 
number of real parameters in B is simply the sum of the 
parameters in each block. It should be clear that H and 
K each contribute L real parameters. The number of real 
parameters in an L x L unitary matrix is so the B-rep 
has + 2L real parameters, the same as the U-rep. We 
said earlier that M-reps require 3L^ -I- L real parameters, 
which is 2L^ — L more than the B-rep implying that (|123p 
is one-to-many from *B to 9Jt. However, observe that for 
any 2L x 2L orthogonal matrix 

M = \/lHSt(yQ, -iVQ)0 (131) 

is not only a valid M-rep but its RHS has + L num- 
ber of real parameters, exactly the number of parame- 
ters required to specify M. It is then very natural to 
ask whether every valid M-rep can be factorized in the 
form of (I13ip . There does not seem to be a simple proof 
in general (for any L) but we do analyze the simplest 
(L = 1) instance of (I13ip in Appendix [Cl In summary 
our analysis shows that for L — 1 and for any M we 
can find a (B,0) such that (|13ip is satisfied. However 
it involves the solutions of transcedental equations. The 
appearance of these transcedental equations, even in the 
simplest case, suggests why proving (|13ip in general may 
be nontrivial. But based on the parameter-count and the 
verification of the L = 1 case we conjecture that every 
M E DJl has a matrix decomposition in the form of (|131l) . 
This motivates us to define 

<Bo = {(B,0)|Be<B,0^ = 0-i} . (132) 

Equation (I13ip is then shown in Fig.[2]by the dotted line. 
This is a many-to-one map from *Bo to DJl as there will 
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be more than one orthogonal matrix that will factorize 
M. This can be seen directly by considering the L — 1 
case which we have analysed in Appendix O From Fig. [7] 
we can see that given an M (and therefore r and 6) there 
is a range of values for ip (which in turn determine 
via (|C1|) ) at one's disposal for factoring M in the form of 

(HSU. 

We can connect the orthogonal matrix in (|13ip and 
(I132p to Fig. m Since given an M e 9Jl we can write 
M' = MO which is also in (|M1) and imply that 

(yM') = (yMo)=0^(yM) = 0^(yB). (133) 

That is to say, if ye emulates y^, then O^yb (obtained 
from post-processing ye) will emulate yivio- 

The completeness of the B-rep also means that every 
valid U must have a matrix decomposition of the form 
given in (in5|) - P^ . The map «B ^ it is thus onto. 
Even though B and U are parameterized by the same 
number of real parameters the map !B — >■ il is actually 
many-to-one. To see this we can consider (I125p - (|128l) for 
L = 1. In this case U is 2 x 2 and S = e*"^, H = 77, and 
Q = 6* as in Sec. IIVB 11 It is sufficient to consider ij — 1, 
in which case (jl25l) - (|128p simplify to 

Uii 61 cos2</>-h (1 - 6*) sin^^ , (134) 
U12 = (1 - 26*) cost/) sin0 , (135) 
U22 sin^ (j)+{l-e) cos^ , (136) 

where U12 is identically equal to U21 as can be verified 
directly from (|126l) and (|127l) for any L. Remembering 
that cos((/)-|-7r) = — cos (j) and sm{(j)+TT) — — sin 0, we can 
see that Bi = (l,e**,6') and B2 = (1, e*('^+^\ 6*) (recall 
the order of the B-rep triple from (|109|) ') both map to 
the same U but yet Bi 7^ B2. We can consider another 
example in which B3 = (1, e*('^+''/^\ 1-6') and Bi map to 
the same U but Bi ^ B3 (this time noting cos{(j) + tt / 2) — 
— smcj) and sin(</) + 'k/2) = cos cj)). The relation between 
03 and il is shown in Fig. [2] 

V. CONCLUDING REMARKS 

A. Summary of Key Results 

In this paper we have put forth two new parameteri- 
zations of diffusive quantum measurements. The first is 
the M-rep introduced in Sec. IIIIl defined by ([7^ . This 
is a characterization of diffusive measurements in the 
Schrodinger picture, and as a consequence is a rather 
abstract parameterization. However it embodies all dif- 
fusive measurements by a single equation which relates 
simply to the efficiency with which each output field is 
measured a/H . The M-rep has been shown to be a com- 
plete parameterization — every diffusive quantum mea- 
surement can always be associated with an M e 931, and 
every M G 9Jt must correspond to a diffusive measure- 
ment. 



We then turned to the Heisenberg picture in Sec. IIVI 
and proposed a realization for arbitrary diffusive mea- 
surements. This gave another parameterization which 
we have called the B-rep (or block-rep) , defined by (|109p . 
By construction every B £ S must correspond to a dif- 
fusive measurement. In particular we derived (jl3ip and 
used it to investigate the existence of B-reps for arbitrary 
diffusive measurements. An analysis of (|131l) for L = 1 
reveals that it requires solutions to transcedental equa- 
tions, which do indeed exist, for any M. The occurrence 
of transcedental equations even in this simplest case sug- 
gests that proving (jl3ip for any L may be non-trivial, or 
even formidable. The verification of (|13ip and the fact 
that this equation (|13ip has the same number of param- 
eters on each side prompts us to propose it as a general 
matrix decomposition for every M e 931, and thereby con- 
jecturing the completeness of the B-rep. 

Finally, we have explained the relations between the 
different parameterizations (including the U-rep) and 
Fig. [2] provides a summary of these relations. 



B. Future Prospects 

Quantum measurements are an integral part of many 
quantum computing architectures. The monitoring of 
qubits (or qunats [59j) and their unwanted environmen- 
tal coupling provides a natural setting to apply quantum- 
trajectory theory. Indeed, the applicability of quantum 
trajectories to quantum information was recognized in 
1995 by Schack and colleagues jl^l- At the same time 
a second group (Pellizzari et al.) had already applied 
quantum trajectories to study decoherence in the cavity 
QED architecture of quantum computing [6l|. Subse- 
quently quantum-trajectory theory was adopted to study 
the effects of decoherence on the quantum information 
processing of discrete-level systems to some generality in 
Ref. 16 3, an d more specifically in quantum error correc- 
tion j63l - (66j , gate performance [67| , and teleportation (of 
both discrete and continuous variables) |68l - [7g |. 

As we scale up our quantum processors one would ex- 
pect quantum trajectories to remain useful. A scale-up of 
existing protocols would inevitably mean a scale-up of the 
necessary measurement from being a few-port device to 
many-port ones. In fact, the error syndrome which condi- 
tions the corrective displacement of just one qunat in the 
nine-qunat error correction scheme of Braunstein (tiI [72j 
already requires eight homodyne detectors. If one wanted 
to simulate the conditioned evolution based on this eight- 
homodyne output using (|24l) . then it is necessary to pa- 
rameterize this measurement as an M , in this case of di- 
mensions of 8 X 16. If all we knew was MM'^/fi e then 
it is not so easy to see what M should be, at least not 
without some thought. Parameterizing this measurement 
into a B-rep on the other hand is much easier. One could 
then use (I13ip to obtain the appropriate M. This illus- 
trates that (jl31l) can be useful when given a (B, 0) € *8o- 
Given that simulating (pure-state) quantum trajectories 
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requires less time and computer memory than a direct 
integration of the corresponding master equation [J] one 
can expect apphcations of quantum trajectories to quan- 
tum information to continue in the foreseeable future. 
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Appendix A: VECTOR-OPERATOR ALGEBRA 

1. Basic Operations 

In the same manner that we write a set of equations 
compactly using vectors and matrices in classical theo- 
ries, we would also like to write a set of operator equa- 
tions compactly. In order to do this we have to introduce 
a vector whose elements are operators which we call a 
vop, short for vector-operator. This appendix will make 
clear particular definitions that we adopt for vops and 
their algebra. 

We begin with the definitions of the vop, its transpose, 
and its Hermitian conjugate. For an arbitrary vop A 
consisting of n components, 



A2 



= Ul,A2,. 



(Al) 



we define its Hermitian conjugate as 



(A2) 



Note the Hermitian conjugate of A is defined by trans- 
posing A (as usually defined), and then taking the Her- 
mitian conjugate of each element of A. We will say that 
a vop is Hermitian if and only if A^ — A^. Matrix- 
operators also appear in this article and we will define 
the Hermitian conjugate of an m x n matrix-operator A 
as follows: 



A = 



/ All A12 
A21 A22 



ml 



At 



/ A^ A^ 
/ '^ii ^21 

At At 



V At 



■ Ai„ \ 

At \ 

At / 



(A3) 



(A4) 



We will denote, for this appendix only, matrix-operators 
(which we refer to as mops), by San-serif letters with a 
hat on top. We will also refer to a 1 x 1 operator as a sop, 
short for scalar-operator. We retain the familiar proper- 
ties of the Hermitian conjugate for complex matrices: 



(A^)t.(At)\ 
(Ab)^ = BtAt , 



(A5) 
(A6) 



where A and B are respectively m x n and n x k mops. 
Equation (|A5[) crops up often so for convenience we define 



A* 



(A7) 



It will also be useful for us to define the real and imag- 
inary parts of an arbitrary sop c to be the Hermitian 
operators 



5R[c]^i(c + ct), 3[c] 



1 

Yi 



(c-ct). 



(A8) 



Since we also deal with complex matrices, we have also 
defined the real and imaginary parts of a complex matrix 
C to be 

5R[C] = i(C + C*), 3[q = i(C-C*), (A9) 
2 Zi 

where C* denotes the complex conjugate of C. We gen- 
eralize the definitions of the real and imaginary parts of 
a sop to a mop: 



3? [A] 



/3?[Aii] 3fi[Ai2] 



mlj 



A-l-A^ 



mm] \ 



/ 3[Aii] 5[Ai2] • • • 5[Ai„] \ 



[A] 



S[A2l] 



\s[A™i] 



[Amn] / 



2i 



(AlO) 



(All) 



(A12) 



(A13) 



Note the necessity of taking both the transpose and Her- 
mitian conjugate in (|A10I) and (|A12I) in order to retain 
an m X n mop. 

Multiplication of A by an arbitrary sop c is defined 
component-wise by 



;A = 



^ cAii cAi2 • • • cAi„ ^ 
CA21 CA22 



\cA 



(A14) 



ml 



7 
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and similarly for the action of an arbitrary superoperator 
AT, 



AfA = 



7VA21 7VA22 



\Afk 



ml 



(A15) 



It is straightforward to see (and will be useful to remem- 
ber), from (|A14|) that 



cA^ = {ckf , 
(CA)^ = A^C^ , 



(A16) 
(A17) 



for any k x m matrix C. 

Two special cases of (IA6|) that we are interested in are 
the inner- and outer-products between vops for which 
[6l) reads 



(A^B)^ = BtA* 
fABT)^ = B*At 



(A18) 
(A19) 



The following identities will also be useful for manipulat- 
ing vop products 



(A^B)C = (A^BC^ 



AB 



A' (BC 



(A20) 
(A21) 



The tensor product for vops are defined analogous to 
standard multiplication rules. If A and B have different 
dimensions then we define 



A«)B 



/ill 
A2 I 



• Bi A2 



B2 
B2 



\Am^Bi 



A^(i)Bn) 



When A and B have the same dimension then 
A^ (g B = ^ ifc (g) Bk . 



(A22) 



(A23) 



We will assign the following symbols for identity oper- 
ations with respect to multiplication: 

• 1 for a sop identity, 

• Ijv for an N X N identity matrix, 

• Ijv for a mop identity (defined by Ijv = Ijvl), 

While it is common to find special symbols for identity 
operations for multiplication in the literature, it is also of 
common practice to not denote the identity with respect 
to addition by any distinguished symbol except for 0. We 
will follow this convention as well. 



2. Commutativity 

a. Sop-brackets 

Unlike normal vectors the order of the vops in an inner- 
product may not be changed in general, 

B^A ^ AJB , (A24) 

and this motivates us to define the sop-bracket 

[A, Bj = A^B - B^ A , (A25) 

so called because it maps two vops to a sop. We then have 
the following sufficient (but not necessary) condition for 
reordering the inner-product. 



Vfc, [Ak,Bk]=0 => [A,Bj =0. 



(A26) 



That [Ak,Bk] = is not necessary for [A,BJ =0 can 
be seen in the following example, 

A^{q,pf, B = (p,g)^, [q,p]=th. (A27) 

Here Ai does not commute with _Bi, and A2 does not 
commute with B2 but we still have A^B = B^A. Note 
that we have not overloaded the usual notation for com- 
mutator brackets in (jA25p because while the sop-bracket 
satisfies the following (which is easy to check) 

[A,BJ = - [B,AJ , (A28) 
[A + B, Cj = [A, Cj + [B, Cj , (A29) 

it cannot satisfy a property analogous to the Jacobi iden- 
tity which the standard commutator does. 

The inner-product is a sop so by definition A^B = 
(A^B)^. Thus (|A26p may also be seen as a sufficient 
condition for (A^B)^ = B^A. 



b. Mop-brackets 

For the vop outer-product we do have a necessary and 
sufficient condition for writing the transpose of AB^ as 
a product of transposes: 



Vi,fc, =0 



AB 



Therefore we define the mop-bracket 
[A,B] = AB^ - (BA^ 



B A^ . (A30) 



(A31) 



which maps vops to mops, and refer to A and B as com- 
muting vops if and only if [A, B] =0. It is easy to see 
that this satisfies 

[A,B] = - [B,A]^ , (A32) 
[A + B, C] = [A, C] + [B, C] . (A33) 
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The symbol for the mop-bracket in (jA31[) has the ceihng 
(] ) on B to remind us of AB^, which is a mop. For the 
sop-bracket of (IA25[) the ceihng is placed on A to remind 
us of A^B. Note the commutativity of vops thus defined 
means that an arbitrary vop will not necessarily commute 
with itself as that is a question about its components. 

From the above one may verify the following formulae 
which help simplify expressions containing mop-brackets: 
For any matrix C (with the appropriate dimensions) 



For any sop c, 



[CA,B] =C[A,B] , 
[A, CB] = [A, B] . 

[c,A] = [cX] , 
[A,c] = [A, a] . 



(A34) 
(A35) 

(A36) 
(A37) 



Note that it is the standard commutator which appears 
on the RHS of (jASGl) and (|A37l) . We also have 



[Ac,B] = A[c,B^] + [A,B]c 
[cA,B] =c[A,B] + ([c,B]A" 

[A,Bc] = (b[aT,c])^ + [A,B]c 
[A,cB] =c[A,B] + [A,c]B^ . 

3. Trace and Average 



(A38) 
(A39) 

(A40) 
(A41) 



Since both matrices and operators appear we will de- 
note an operator-trace by Tr[c] and the trace of a matrix 
C by tr[C]. This distinction is especially important for 
operators with a continuous eigenvalue spectrum where 
the operator-matrix analogy breaks down. The trace of 
an m X n mop is a matrix in C™^" 



Tr[A] 



/Tr[An] Tr[Ai2] 
Tr[A2i] Tr[A22] 



Tr[Ai„] \ 



V Tr[A™i] 

From this it should be simple to see that 

Tr[AT^(Tr[A])\ 
(Tr[A])* =Tr[A*] . 



(A42) 



(A43) 
(A44) 



Particular cases of (|A43p and (IA44I) which we are inter- 
ested in are 



:Tr[A] 



Tr[At] , 



(Tr[ABT])* =Tr[A*Bt] 



(A45) 
(A46) 



If A is any m x n mop and B is n x fc, then we have the 
"cyclic" property 



Tr[(AB)"^] = Tr[B^A^] 



(A47) 



Special cases of the "cyclic" property (|A47I) are 

Tr[cA] = Tr[Ac] , (A48) 
Tr[(AB^)^] = Tr[BA^] . (A49) 

On some occasions we will want to permute a product 
of three vops inside a trace and the following identity 
may be verified, 

Tr[AB^C] =Tr[B^CA] = Tr[(CA^)^B] . (A50) 

An important tool from quantum statistics that we 
will use is the quantum regression formula, also known 
as the quantum regression theorem, formulated by Lax 
[zIjItI]. a good discussion of this can be found in [sJl- In 
essence this result expresses the average of a product of 
two operators, each at a different time, as a trace in the 
Schrodinger picture. To write down the vop counterpart 
of this result we first recall the regression formulae for 
sops: Given any two operators A and B in the system 
Hilbert space, and the solution to the Markovian master 
equation p = Cp, 

{A{t) B(t + t)) = Tr{B(0) e^" [p(t) i(0)] } , (A51) 
{A{t + t) B{t)) = Tr{i(0) e^" [b(0) p{t)\ } . (A52) 

for any r > . Let us define the quantum average of a 
mop A at time t as 

(A) - Tr[p(i)A] . (A53) 
A corollary which follows from this definition and (|A45p 



IS 



(cA)t = (Atct) 



(A54) 



obtained by letting A — > pcA in (jA45p . We can now 
extend (|XgT|) and ((X52)) to 



{k{t) B^(i + r)) = (Tr{B(0) e^^ [p{t) A^(0)] }) 



(A(i + r) B^{t)) = Tr{A(0) e^" [b^(0) p{t)\ } 



(A55) 



(A56) 



4. Superoperators of Special Interest 

When considering continuous measurements terms of 
certain forms arise frequently. Here we will define super- 
operators whose forms allow us to write these frequently 
occuring terms compactly. For an arbitrary sop s we 
define 



J\A]B = ABA^ 



(A57) 



Note that A here should be treated as a parameter and 
B the input for i7[j4], the output of the superoperator is 
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another sop given by the RHS of (IA57[) . The definition 
(|II7)) thus impUes that J[A^B = A'f B A. 

Often we encounter a sum of terms in the form of ()A57P 
so we define 

J[k]B = J[Ak]B ^Y.AuBAl= k^B A* . 

fc k 

(A58) 

From (|JT5t and (IX58| the action of J' [A] on a vop pro- 
duces another vop, given by 

:7[A]B = ^ifcBit =(AB^)^A* . (A59) 

k 

This also gives J'[A-'-]B — (A^B^)^A . Here we note 

that J7[A-'"] is formally equivalent to its superoperator 
adjoint [ sH- 

From (jA57p we can define what is sometimes referred 
to as the "dissipator" [sl] 

V[A]i3 = J[A]B-]^{A'^A,B], (A60) 

where {A, B] = AB + BA. When A ^ A^ one may 
prefer to write the dissipator as a nested commutator, 

V[A]B^~\[AM:B]]- (A61) 

As in (|A58p . we shorthand a sum of dissipators by 

V[k]]3 = ^ V[Ak\B . (A62) 
fc 

Equations (|A15|) and (|A62p then permit us to write 

I?[A]B=^[A]B-i{AtA,B} 

= (AB^)^At - i {A^A, B} . (A63) 

In the case when A ~ k^ the vop-generalization of (|A61[) 
is 

I?[A]B--i{(AT[A,B])^+ [B,A]a} . (A64) 

The appearance of a dissipator 2? [A] in the master 
equation is associated with the coupling of the system to 
the environment via A. If the system is also under con- 
tinuous observation then measurement back-action terms 
arise and they can be written concisely by defining 

■H[i]i3 = ii3 + M^-Tr(iS + M^)i3. (A65) 

We will invariably be considering H[C^A] with C = C*. 
In this case one may prefer to generalize (|A65P to 

'H[k]B = kB + Bk^ - Tr[AB -I- BA^] B (A66) 



and write 

■H[C^A] = C^H[A] . (A67) 

Note that (|A65p and (|A66[) each contain a term which is 
nonlinear in B. One may find it useful to also define the 
linear versions of (|A65|) and (|A66p . 

1l[A]B = AB + BA'^, (A68) 
'H[k]B = kB + Bk^ . (A69) 

Appendix B: DERIVATION OF Eq. (|46|l 

We first write (suppressing the time-dependence) 

Tr[(p-f dpe)'] ^l+'^Ty[p{hdp,)\+^Tv[{hdp,f] , 

(Bl) 

where the change in the state is given by 

ftdpc + 'H[dw^M^c]p . (B2) 

Recall that C is defined by ([1]). We will assume that 
p{t) is unconditioned and such that Tr[p(i)] = 1 and 
p\t)^p{t). 

We first work out the second term in (|Bip . Taking the 
trace of (|B2p against p, 

Tv{p{hdpc)] Tr{p/:p}rfi + Tr{pH[rfw^Ml'c]p} . 

(B3) 

From (|M5)) we find, for any A, 

i:v{p-H[A]p] =Q . (B4) 

From ((!]) we get 

Tr{p/:p} = (ct)(c)-(ctc) . (B5) 

Therefore (|B3| is 

Tv{p{hdpc)] = Tv{pCp]dt = ((c^)(c) - (c^c))dt . 

(B6) 

To order dt the third term in (|Bip is proportional to 
Tv[{hdp^f] = Tr[(-H[dw^M^c]p)^] . (B7) 
For any A, 
{H[A]pf 

- {Ap){Ap) + {Ap){pA^) - {Ap)Tv[Ap + pA^]p 
+ + {pA^){pA^) - {pA^)Tv[Ap + pA^]p 

+ Tr[ip + pi^]Tr[ip + pit]p2 . (b8) 

Taking the trace and using p^ = p, 

Tr{(H[i]p)'} = 2((iti) - . (B9) 
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This gives, for A = dw^M^c, 

Tr{(-H[dw^Mtc]p)^} = 2((c^MMtc) - (c) MM\c)) dt . 

(BIO) 



Thus substituting dBe]) and (IbTO)) into (|BT]) . 
Tr[(p + dpe)'] 

= 1 + ^tr[((c- (c))*(c - (c))^)(l, - MMyhf]dt 
= 1 + tr[H {l^ - MMy ]dt . (Bll) 
The matrix H was defined in (|T7)) . 

Appendix C: Eq. p3T|) FOR L = 1 

We can show that for L = 1 it is always possible to find 
a decomposition in the form of (jl3ip for an arbitrary 
M. The most general 2x2 orthogonal matrix can be 
parameterized by a single variable ip, 



= 



cos ip sm ip 
=p sin If ± cos if 



(CI) 



where the sign in the second row corresponds to det(O) 
±1. This gives, from ([m]) . 



If we define 



cos if sm ip 

sm ip 

cos p 



mi 
m2 



\mi\e 
|r7i2|e* 



then our problem is to find (77, (/), 9, p) such that 

|TOi|e*"i = V^e*'^(^ cos p±iV^ sin(p) 
|rn2| e'"^ = \/hrje^'^{V6 simp q= cos 1^9) , 
given mi, m2, which determine 77 by 



\mi\ 



(C2) 

(C3) 

(C4) 
(C5) 

(C6) 



Equations (IC4[) and (IC5[) are true if and only if the mod- 
ulus and argument (phase) of each side are equal, i.e. 



|mi p = hrj [9 cos^ p + 9 sir? ip] 
|m2p = h'q[9siv? ip + ^cos^ (p] 



and 



ai = (j) + arg [\f9 cos ip ± i \/^ sin 93) 
a2 = + arg (v^ sin =F cos 
Rearranging (jC7|) we find 

|toiP — /iTysin^ 



fi?/ cos(2<y9) 



(C7) 
(C8) 

(C9) 
(CIO) 

(Cll) 



If we now substitute (jC6|) into (jCll|) and define r 
(|mi|/|m2|)^, then we get 



.{lp) = 



2 • 2 

r cos ip — sm 

(r + l)(cos2 ip — sir? ip) 



(C12) 



Note that (p must also ensure < < 1. For small r we 
find, from (ICT2l) . 



9o('^) 



• 2 

■ sm ip 



9 • 2 
cos-^ (p — sm 



while for r large 9r{y^) approaches 
9oo{v) 



COS^ iy9 



cos^ f — sin^ </? 



(C13) 



(CM) 



At these extremes the only values of ip that enforce 9 G 
[0, 1] are integer multiples of 7r/2, for which 9 is either 
or 1. 

If we now attempt to eliminate 9 in (jC8[) by substitut- 
ing in (jC12p we simply arrive at (jC6|) . which is indepen- 
dent of ip. Thus (jC7p and (jC8[) are solved simultaneously 
with any ip and (jC12p . The value of ip is determined from 
solving (IC9p and (|C10p simultaneously. First eliminating 
4> we see that (/3 must satisfy 

S = ai-a2 = F^{^) , (C15) 

where we have defined the argument of mi/m2 as 



F^{ip) = arg 



\j9r{'p) cos 93 ± i\/l — 9r{ip) simp 
^ 9r{(p) sin (p ^ Z^l — 9r{ip) COS ip 



(C16) 



We have written 9r{'p) in (jC16P to signify that here we are 
substituting in (IC12p . Recall the signs here correspond 
to det(O) — ±1. Equation (jClSp is transcedental. We 
plot F^{ip) in Fig. [5] as a function of r and ip. Wc find 
that for det(O) = 1 



V9 



cos ip 



sm ip 



\f9 sin p — — 9 cos ip 
9 sin^ ip + {1 — 9) cos^ p 



> 



(C17) 



which means the range of should be between and tt 
as seen in Fig. |6l A noteworthy feature occurs at r = 1. 
From (IC12I) we get 9i = 1/2 which gives 



^1^ 



arg 



cos ip + I sm p 
sin ip — i cos ip 



arg[i] ^ 7r/2 . (CIS) 



We also plot F:^[ip) (sohd blue line) together with 
9r{ip) (purple dashed line) for selected values of r as a 
function of ip in Fig. [71 Observe that F^{ip) always has 
a support which corresponds to 6*^ € [0, 1]. We illustrate 
this region in Fig. EKb) with a light shade. It can be seen 




FIG. 7; Plots of F^{(p) (solid blue line) and 9r{(p) (dashed purple line) for selected values of r: (a) r = 0.05, (b) r = 0.5 (c) 
r = 2, (d) r = 20. 



21 



in Fig. [7] that this region persists for the various vahies 
of r that we have chosen. This ihustrates that there is 
always a value of tp which solves (jClSP for 5 £ [0, tt], and 
for which <E [0, 1]. Note that more than one value of tp 
will solve (jClSP for a given 5 which means that M can 
be factorized by more than one (recall (jCip and (|C2p ). 
Having obtained a value of ip we can then substitute it 
back into either (|C9|) or (jC10|) to obtain (j). 

In the above we have concentrated on the case 6 e 
[0,7r]. Here we show how M can be factorized by (|13ip 
if we were given S G [— 7r,0]. The analysis up to (|C15|) 
would remain the same but in this case we can first let 
a2 = Q^2 i and solve 



where 



|mi|e" 



\mi\e" 



= y hij e 



is factorized by 



,0 ( costs' - simp' 
sin p' cos p' 




Substituting this into (IC20I) and taking the transpose we 
see that M has the form of (|13ip with the orthogonal 
matrix given by 



ai 



(C19) 



for p' . The factorization of is then given by 







cos p sm p 
— sin p' cos p' 



1 

-1 



(C22) 



|mi|e" 
I male'"' 



|m2|e*KT^) 

1 \ /|mi|e*"i 
-1 j Vl"^2|e'"^ 



where p' defined by (IC19I) . Note that (IC22I) now has a 
determinant of minus one. We see then for 5 G [— tt, 0] 
(C20) M is stih factorized by but with an such that 

det(O) = -1. 
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